RIGOROUS POINTWISE APPROXIMATIONS FOR INVARIANT 
DENSITIES OF NONUNIFORMLY EXPANDING MAPS 
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Abstract. We use an Ulam-type discretization scheme to provide pointwise 
approximations for invariant densities of interval maps with a neutral fixed 
point. We prove that the approximate invariant density converges pointwise 
to the true density at a rate C* ■ , where C* is a computable fixed constant 
and m~ 1 is the mesh size of the discretization. 



1. INTRODUCTION 

Ulam-type discretization schemes provide rigorous approximations for dynamical 
invariants. Moreover, such discretizations are easily implementable on a computer. 
In [TB] it was shown that the original Ulam method [33J is remarkably successful 
in approximating isolated spectrum of transfer operators associated with piecewise 
expanding maps of the interval. In particular, it was shown that this method 
provides rigorous approximations in the L 1 — norm for invariant densities of Lasota- 
Yorke maps (see [18] and references therein). This method has been also successful 
when dealing with multi-dimensional piecewise expanding maps [20] . and partially 
successful Q in providing rigorous approximations for certain uniformly hyperbolic 
systems jTU] HJ. Recently, Blank [5] and Murray [2T] independently succeeded 
in applying the pure Ulam method in a non-uniformly hyperbolic setting. They 
obtained approximations in the L 1 — norm for invariant densities of certain non- 
uniformly expanding maps of the intervaj^j 

Although L 1 approximations provide significant information about the long-term 
statistics of the underlying system, they are not helpful when dealing with rare 
events in dynamical systems. In fact, when studying rare events in dynamical 
systems [TJ 115] one often obtains probabilistic laws that depend on pointwise in- 
formation from the invariant density of the system. In particular, extreme value 
laws of interval maps with a neutral fixed point depend pointwise on the invariant 
density of the map [13] , 

Statistical properties of non-uniformly expanding maps were studied by Piani- 
giani in [22] who first proved existence of invariant densities of such maps. Later, it 
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^See [6] for examples where the pure Ulam method provides fake spectra for certain hyperbolic 
systems. 

2 In 1211 . in addition to proving convergence, Murray also obtained an upper bound on the rate 
of convergence. 
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was independently proved in [14] [19j [24] that such maps exhibit polynomial decay 
of correlations. The slow mixing behaviour that such maps exhibit has made them 
good testing tools for real and difficult physical problems. 

The difficulties in obtaining pointwise approximations for invariant densities of 
interval maps with a neutral fixed point is two fold. Firstly, the transfer opera- 
tor associated with such maps does not have a spectral gap in a classical Banach 
space. Therefore, powerful perturbation results |16j^| are not directly available in 
this setting. Secondly, invariant densities of such maps are not L°° functions. Con- 
sequently, to provide pointwise approximation of such densities, one should first 
measure the approximations in a 'properly weighted' L°°— norm. 

In this note we use a piecewise linear Ulam-type discretization scheme to provide 
pointwise approximations for invariant densities of nonuniformly expanding interval 
maps. We prove that the approximate invariant density converges pointwise to the 
true density at a rate C* ■ where C* is a computable fixed constant and to -1 
is the mesh size of the discretization. To overcome the spectral difficulties and the 
unboundedness of the densities which we discussed above, we first induce the map 
and obtain a uniformly piecewise, expanding and onto map. Then we perform our 
discretization on the induced space. After that we pull back, both the invariant 
density and the approximate one to the full space and measure their difference in 



a weighted L°°-norm. Full details of our strategy is given in subsection 3.2 



In section [2j we recall results on uniformly piecewise expanding and onto maps. 
Moreover, we introduce our discretization scheme and recall results about uniform 
approximations for invariant densities of uniformly piecewise expanding and onto 
maps. In section [3j we introduce our non-uniformly expanding system, set up our 
strategy, and state our main results, Theorem 3.1 and Corollary |3.2| Section [4] 



contains technical Lemmas and the proof of Theorem |3.f 

2. Preliminaries 

2.1. A piecewise expanding system. Let (A, 05, A) denote the measure space 
where A is an interval, *B is Borel u-algebra and A is normalized Lebesuge measure 
on A. Let T : A — > A be a measurable transformation. We assume that there exists 
a countable partition V of A, which consists of a sequence of intervals, V = {Ii}^ , 
such that 

(1) for each i = 1, . . . , oo, T; L :— To is monotone, C 2 and it extends to a C 2 

function on If, 

(2) Ti(Ii) = A; i.e., for each i = 1 , . . . , oo, Ti is onto; 

(3) there exists a constant D > such that sup, sup xeI . jf,^[ — D i 



(4) there exits a number 7 such that jiy < 7 < 1. 



Let £ : L 1 — > L 1 denote the transfer operator (Perron-Frobenius) [4] [8j associated 
to f: 

^See also 1 1 21 for another perturbation result, which also requires a spectral gap. 
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Under the above assumptions, among other ergodic properties, it is well known (see 
for instance [7]) T admits a unique invariant density /; i.e. Cf = /. Moreover, C 
admits a spectral gap when acting on the space of Lipschitz continuous functions 
over A We will denote by BV(A) the space of functions of bounded variation 
defined on the interval A. Set || • ||bv(A) : — Va + || ■ Hi, where Va denotes the 
one-dimensional variation over A. Then it is well known that (BV(A), \ \ ■ \\bv(A)) 
is a Banach space and C satisfies the following inequality (see [22] for instance): 
there exists a constant C^y > such that for any / G BV(A), we have 

(2.1) V A Cf < jV A f + C LY 



Inequality (2.1) is called the Lasota-Yorkc inequality. 



2.2. Markov Discretization. We now introduce a discretization scheme which 
enables us to obtain rigorous uniform approximation of / the invariant density of 
T. We use a piecewise linear approximations which was introduced by Ding and Li 
9J. Let rj = {ci}"L be a partition of I into intervals. Since uniform partitions are 
the first choice for numerical work, we set c, — Cj_i = Everything we do can be 
easily modified for non-uniform partitions with only minor notational changes. Let 



<Pi = X[ Ci _i, Ci ] and <j>i{x) = ml ipidX. 

Jo 

Let ipi denote a set of hat functions over 77: 

(2.2) Vo := (1 - <t>i) Am ■= 4>m and for i = 1, . . . , m - 1 , ipi := (fa - cj> i+ i) 
For / e L 1 , we set Jj := [cj_i,c,-] and 

fi :=m J f dx, i = 1,2, . . . m, 
the average of / over the associated partition cell. For / G L 1 we set 



m — 1 



if ■= .Ma + ^2 ~ ^ + 



i=l 

Obviously, the operator Q m retains good stochastic properties; i.e., 

• for / > 0, Q m f > 0; 

• jQ m f = Jf- 

We now define a piecewise linear Markov discretization of C by 
(2.3) P m := Q m o C. 

Notice that P m is a finite-rank Markov operator whose range is contained in the 
space of continuous, piecewise linear functions with respect to rj. The matrix rep- 
resentation of P m restricted to this finite-dimensional space and with respect to the 
basis {ipi} is a (row) stochastic matrix, with entries 



Pij := m / Cipi > 0. 



In [2], a Lasota-Yorke inequality was obtained for Markov interval maps with a finite partition. 
The proof carries over for pie cewise onto maps with a countable number of branches satisfying 
assumptions of subsection 2.1 
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By the Perron-Frobenius Theorem for stochastic matrices [17] , P m has a left invari- 
ant density / m ; i.e., 

fm fm^m • 

The following theorem was proved in [2]: 

Theorem 2.1. There exits a computable constant C such that for any to G N 

/ — fm \\oo 5: C . 

TO 

Remark 2.2. We recall that in [5] it was shown that the constant C, which is 
independent of to, can be computed explicitly. 

3. Pointwise Approximations for invariant densities of Maps with a 

neutral fixed point 

3.1. The non- uniformly expanding system. Let I = [0, 1] be the unit interval, 
A be Lebesgue measure on [0, 1]. Let T : I —> I be a piecewise smooth map. We 
assume that 

• T(0) = and there is a x G (0, 1) such that T x = T | [0)Xo) , T 2 = T \ [xoA] 
and T x : [0,x ) [0, 1),T 2 : [x , 1] °4° [0, 1]; 

• Tx is C 1 on [0,x ], Ti is C 2 on (0,x ] and T 2 is C 2 on [x , 1]. 

• T'(0) = 1 and T'(x) > 1 for x G (0,x ) ; \T'(x)\ > /3 > 1 for x G (x ,l); 

• Ti and T[ have the form 

Ti(a:) = s + x 1+a + x 1+a 6 (x), 

T[(x) = 1 + (1 + a)a; Q + x a S 1 (x), 
where, < a < 1 and (5j(x) — > as x — > for i = 0, 1 with > 0. 

It is well known that T admits a unique invariant density /* []3J [THl HU [Mj and 
the system (7, T, /* - A) exhibits a polynomial mixing rate [HI HH1 121] ■ Moreover, it 
is well known [TH[Tni[21] ^ na ^ ^ ne T-invariant density, /*, is not an T°°-function. 
In particular, near x = 0, f*(x) behaves like x~ a . Despite this difficulty, we will 
show that, for any x G (0, 1], one can obtain rigorous pointwise approximation of 
/*(*)• 

3.2. Strategy and the statement of the main result. We first define a Banach 
space which is subset of L , but where /* has a finite norm. More precisely, let B 
denote the set of continuous and integrable functions on (0, 1] with the norm 

||/|| B = sup \x 1+a f(x)\. 

ag(0,l] 

When equipped with the norm || • B is a Banach space. The fact that /* G B 
follows from Lemma 3.3 of [14] . Our strategy for obtaining pointwise approximation 
/* consists of the following steps: 

(1) We first induce T on A C I and obtain a T which satisfies the assumptions 
of subsection 12. II 
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(2) On A, we use Theorem 2.1 to say that f m , the invariant density of the 
discretized operator P m := Q m o C, defined in equation (2.3), provides a 
uniform approximation of / the T-invariant density. 

(3) Then we write /* in terms of /, and define a function f m as the 'pullback' 
of f m - 

(4) We then use steps (2) and (3) to prove that ||/* - f m \\ B < C*^, and 
deduce a pointwise approximation of /* . 

3.2.1. The induced system. We induce T on A := [xq, 1]. For n > we define 



Set 



W := (x ,l),and W„ ■= (x n ,x n -i), n>l. 
For n > 1, we define 

Z n := T 2 "" 1 (W„_i). 
Then we define the induced map T : A — > A by 
(3.1) f (x) = T n {x) for x e Z n . 

Observe that 

T(Z n ) = W n -i and r Zn = n, 
where Tz n is the first return time of Z n to A. An example of the map T and its 
induced counterpart T are shown in Figures flland [2] respectively. It is well known 
(see for instance [53]) that the T defined in ( |3.1[ ) satisfies the assumptions of sub- 
section [2~l] and, by Theorem |2.1[ one can obtain a rigorous uniform approximation 
of its invariant density /. Moreover, by Lemma 3.3 of [3], /*, the invariant density 
of T, can be written in terms of /: 




0.2 0.4 0.6 0.8 1 



Figure 1. A typical example of a map T which belongs to the 
family defined in subsection |3.1| 
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Figure 2. This figure shows the induced map t corresponding to 
the map T of Figure [l] 



(3.2) 



c T f(x) 



x) 



for a; € A 



for x e I \ A 



£i Vl Z)T< " ) (T 2 - 1 T 1 - ( ' l - 1) I 
where / is the T-invariant density, c^T 1 = Efcli T z k i^{Zk) 1 and fi = f ■ X. 
3.2.2. TTie approximate density and the statement of the main result. Set 

c T , m f m (x) for X e A 

^ I |£)T(")(T~ 1 T 1 _( 'x)| y X 



(3.3) / m (x) : 



def 



where f m = F m f m , and P m is the Markov discretization of £ dehned in (2.3|, 
c r,m = J2T=i T z k fj-rn( z k), and /t m = f m ■ X. We will show that the function / m 
defined in (3.3 1 provides a rigorous pointwise approximation of /*. 

Theorem 3.1. For any m E N we fttwe 

nr-/mii*<c^, 

771 

w/7ere 



C* 



CI 



+ M(l + a) C 4 ; 



in particular, C is the computable constant of Theorem 



2.1 



M :-- 



C := a(1 2 +a) [l + 25 (x ) + 5 2 (x )\, C 1 := (2[2* - I]) 1 /*, 
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and 

C 2 = 1^2^-1]^. 
x 

As a direct consequence of the Theorem [XT] we obtain a pointwise approximation 
of/*: 

Corollary 3.2. For any x <G (0, 1] we have 

,,,, , .... C* In m 

Proof. For a; € (0, 1], we have 

\r(x)-ux)\ = ^i^(/*(x)-/ m (^))i < _l_n r _ /m || B < ^c* 11 ^. 

□ 

4. Proofs 

4.1. Technical lemmas. We first introduce notation of certain functions which 
appear in the proof of Theorem 3.1 For x € I \ A set: 



T>(x) ' 



G\{x) :— — r — ; and for n > 2, G n (x) 



x 1 + a 



Lemma 4.1. For x € I\A, we have 

[1 + x a + x a 6 {x)] 1+a < 1 + (1 + a)[x a + x a S {x)} + " (1 + - [x a + x a 8 a {x)] 2 . 
Proof. Let 

4n(x) := [l + x a +x a S (x)] 1+a 

and 

0a(x) := 1 + (1 + + x^ (a;)] + ^1±^1 [ x « + ^(x)] 2 . 

Note that <f>i(0) — 02 (0) = 1. Therefore, to prove the lemma, it is enough to prove 
that 4>[(x) < <j>2(x). We have: 

^(x) = (l + a)(l + WW 

0' 2 (z) = (l + a)(l + a£(aOK'(*), 

where £(x) := x a + x a 6 (x) > cQ Notice that f'( x) > 0. Thus, we only need to 
show that 

(4.1) (l + £0r)) tt <(l + <(x)). 

Indeed, (|4Tj holds because (1 + £(0)) a = (1 + a£(0)) = 1 and 

[(1 + t(x))<*]' = {1 + ^ x))1 _J '(x) < ai'{x) = [1 + «£(*)]'. 

□ 

5 It is obvious that £(0) = and for x > 0, > 0. 



8 



Pointwisc Approximations of Invariant Densities 



Lemma 4.2. For x E I \ A, we have g(x) < 1 + Cox 2a , where 
Co= a{1 + a) [l + 2So(x )+S 2 (x )}. 



Proof. Using Lemma |4.1[ we have: 

(2jL*)l+a [\ + x a +X a S (x)] 1+a 



g( x ) 



< 



T{(x) 1 + (1 + a)x a +x a 5 1 (x) 
1 + (1 + a)[x a + x a 5 (x)} + ^±^[x a + x a S (x)} 2 

1 + (1 + a)x a + x a S!{x) 
1 + (1 + a)[x a + x a S (x)} ^±[x a + x a S {x)} 2 



1 + (1 + a)x a + x a S 1 (x) 1 + (1 + a)x a + x a S 1 (x) 
< l + ^f^K+^ (x)] 2 

= 1 + a(1 ^ a) (1 + 2S (x) + 5 2 {x))x 2a < 1 + C x 2a . 

□ 

Lemma 4.3. Letx n = T^ n x . Forn> l,x n < C>~s, where C\ = (2[2i - 1]) 1 /". 

Proof. Observe that C\ > 1 > T-, - (o;o) = X\. Therefore, the lemma is true for 
n = 1. Next, for n > 2, we suppose that CE n _i < C\{n — = , and prove that 
%n < Cin~a. If it is false, that is x n > Cin _ °, then by our inductive statement 
on we have: 

Ci(n-l) - « >i„-i =Ti(a; n ) > Cm - - [1 + Cfn" 1 + Cfn^SoiCm-*)]. 
This is equivalent to 

n[(l + -^r)- - 1] > Cf[l + *o(C in -s)]. 

71 — 1 

By convexity of the function z° , it follows -^[2 a — 1] > Cf [1 + <5o(Cin~ ° )], that 
is 

C? < -Ar[2* -l]/[l + <y (Cm-5)] <2[2^ -1] = Cf. 

71—1 

A contradiction. Therefore, a;„ < C\n~~^ 1 and this completes the proof of the 
lemma. □ 

Lemma 4.4. For x G I \ A, we Ziawe 



13 ' 
and /or n > 2, 

G n (x)<M(n-l)-( 1+ i\ 

where M = — — ^ — . 
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Proof. For n = 1, it is easy to see that 



„l+a 



Gi(x) < ^ 



For n > 2, we have 
G n (x) 



|IXr('0(2 , J r 1 ir (B - 1) a:)| 



(4.2) 



^i C^i x ) ' T{ {Ti x) ■ ■ ■ 

( x M+ct ( T i ±x M+a 



-lr r -(n-l) 



("-!)„ 



r{(rr (n - 1) .x) 1^(25-^- 



<g(Tr L x)-g(T^x)---g(T; 



-(«-i) x U+c 







By Lemmas 4.2 and 4.3 for any k > l,x g [0, xo), we have 

S (Tf *(*)) < 1 + CoCir*^)) 20 < 1 + C (T^ k (x )) 2a 
= 1 + C (x k ) 2a <l + C Cf a k- 2 . 



(4.3) 



Therefore, using 14. 2h and (14. 3h , for n > 2, we obtain 



n-l 



G„(x) = l[g(Tr k (x)) 



k=l 
n-l 



{T; {n - 1] {x)) 1+a 

is 



C 1 1+0 (n-l)-( 1+ ^ 



fc=i 



exp{^ln(l + C C 2a fc- 2 )} 



fe=i 

n-l 



Ci +a (rc-l)-( 1+ n) 



< exp{^C C 2Q fc- 2 } 



fc=i 



C 1 1+a (n-l)-( 1 +^) 



n — 1 p 



< M(n- 



Lemma 4.5. 



where C3 = -| 



C 2 



1 /3(l-x ) 



5>-A(Z„)<C 3 , 

n=l 

(a + f=a) and C 2 = ij^+ipj* - 1]*+*. 



□ 



10 



Pointwisc Approximations of Invariant Densities 



Proof. By Lemma 



4.3 



we have X(W n ) = x n _i - x„ = Tx(x n ) -i„ = - 



^Cj +0, n-( 1+ -) = C 2 n-( 1+ H Since T 2 {Z n ) = W n -i, we have 



£VA(Z n ) < + 



n=l 



n=2 

oo 



A(W n _! 

/8 



< 



1 - x n(x w _ 2 - x„_i) 

1 - x y> (n + l)(x n _i - Xn) 

^ n=l K 

1 — Xo , n {Xn-l — X n ) ^ (x n _i — X n ) 

n=l ^ n=l 

oo 



P 



P 



P 



< 



1 - X , C 2 _ 



ft ^ 



c. 



/3 



< 



^ + fd + /^) + f(i + /^ 

i i 

l-x ,C 2 . 2 -a. 

= — 75— + 7r( a + i ) = (1 - so) • C 3 . 

p p 1 — a 

This completes the proof of the lemma since A(-) = j _ x ■ 

Lemma 4.6. FFe /iawe 



Crl <C 3 -C- 



In m 



Proof. Using the fact that c T < 1, c T . m < 1 and Theorem 2.1 we have 



Crl < 



E T Zh p, m {z h ) e T Zk K z k) 

k=l fe=l 
00 

E k[K Z k) - frm(Zk)} 
k=l 



E T"Z fe Am(^fc) ' E T Z k K Z k) 
k=l k=l 



< 



||/-/m||oo (5>*( Z *) 



\fc=l 



< c— -c 3 . 

m 



In the last estimate, we have used Lemma |4.5| 

We now have all our tools ready to prove Theorem |3.1 
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Proof, (of Theorem 3.1) Using (3.2) and (3.3), we have 
(4.4) 

\\f*-fm\\B= SUp \x 1+a (f*(x)-f m (x))\ 

xe(o,i] 

< SUp \x 1+a (r(x)-f m (x))\ + S Up\x 1+a (r(x)-f m (x))\ 
x£l\A i£A 

°o l+ a 

+ sup |a; 1+a (c T /(a;) - c r , m / m (a;))|. 

Notice that for a;6/\A, and n > 1, z„ := T^ -1 !^ ^ x G A. Then using the fact 
that c T < 1, Cr.m < 1, Theorem 2.1 Lemma 4.6 and (4.4), we obtain: 

(4.5) 



„1 + Q 



||/*- /mils < SUp |V X _ f „_ n I' SUp |(c r /(z ra ) - Cr, m fm(z n ) 

xei\A ^-{DT( n )(T^ 1 T 1 ln 1} x) z n eA 



+ SUp \c T f(x) - C T . m f m (x)\ 
x£A 

< sup | ^ 



sup |/(z„) - / m (z„)| + |c T - c T , m | sup |/(z„)l 

z„GA z„SA 



+ sup |/(a;) - f m (x)\ + \c T - c nm | sup |/(a;)| 

igA i£A 



< c 



— I sup V|G„(x)|(l + C 3 sup |/> n )|) + (l + C 3 sup|/(x)|) ] . 

Since / G BV(A), we have sup 2 , eA |/(x)| < U A / + i^||/||i,a- Therefore, using 
the Lasota-Yorke inequality (2.1), we obtain sup^A < + jz^- Using 

Lemma 4.4 and (4.5), we obtain: 



||/* -/mils < C A C— (i+^I + ^MCn-l)-^ 



.1+a 



C 4 C- 



rn 

..lnm 
rn 



„l+a 



n=2 



< CaC 1 



„l+ct 



/3 



M(l 



-(i+s) 



lnm 



□ 
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